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Abstract 



O ' The effective motion equation that describes the different monomer alternation along the heteropolymer 

chain is proposed. On its basis the supersymmetry field scheme that allows to obtain the equations for the 
structure factor and Green function is built up. The memory and ergodicity breaking effects are investigated 
depending on the temperature and quenched disorder of the monomer alternation. The phase diagram that 
, determines the existence of the non-ergodic and freezing states is provided. 

^ ■ PACS numbers: 61.4f .+e, 05.40+j, fl.30.Pb 



1 Introduction 



■ O ■ In recent years investigation of the unusual behaviour of random heteropolymers with a random sequence of 
^ ' different kinds of monomers has attracted considerable attention (see and references therein). With 

temperature decreasing such systems experience two types of phase transitions at which polymer either goes to 
freezing state or has a microphase separation that is inherent in protein biomolecules type. 

At folding process heteropolymer chooses out of infinite set of possible states the unique frozen non- 
equilibrium but stationary state with fixed spatial conformation and monomer sequence (type of DNA). The 
QQ problem in considering of such transitions is to build up more simple method of averaging thermodynamic values 
^\ ' over the quenched disorder. Here, the methods of statistical physics first used in spin glass investigation are 
. taken as basic ones [||. 

Certain of polymers type of block copolymers, undergo the so-called microphase separation at which blocks 
: of different chemical composition segregate to microdomain structures to gain a spatially inhomogeneous com- 
' position throughout the pattern. The chemical bonds between the blocks prevent the usual macroscopic phase 
"Y"; ' separation observed in homopolymer blends. A variety of microphases, such as one-dimensional lamellar struc- 
' tures, hexagonal arrays of cylinders, body-centered-cubic arrays of spheres and so on have been observed for 
^ , simple diblock copolymer melts. 

I At first, microphase separation had been studied within mean- field approximation |^ in the case of A-B 

^ t block copolymers with arbitrary composition /. It turns out that at / 7^ 0.5 with temperature decrease the 
Q ■ system undergoes a sequence of the microphase separations of the first order, where period 211 /k^ of space 
Q , structure is of order of block length and does not depend on temperature. In ^ it was shown a principal role 
^ ■ of the fluctuations in phase transition. By virtue of condition ko ^ the order parameter fluctuations give 
• ^ , divergent contribution to thermodynamic values, so that continuous phase transition transforms to the first 
' order weak transition . 

At passage to random heteropolymers both the microphase separation and the freezing remain, but here 
the space period acquires strong dependence on temperature ||. According to field considerations ^ 
the fluctuations suppress both mentioned transitions in random copolymers, but copolymer melts in disordered 
media can suffer these transformations. 

Similar to spin glass , the main approach of statistical mechanics of random heteropolymers is the replica 
method that is not only pragmatic, but very non-descriptive and artificial. Apart from the replica trick, 
the transfer matrix method kinetic approach and others (see are used in random heteropolymers 
theory as well. However, it is a well known in theory of spin glass that within the framework of Sherrington- 
Kirkpatrick model (l4j both replica method and supersymmetry (SUSY) field approach give identical results 
||l5|| . On the other hand, the advantage of the SUSY field usage is that its components have explicit physical 
meaning |p^. 
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By recent moment, SUSY method in polymer theory had been proposed by Vilgis [|^, but did not obtain 
further development. In this work we intend to complete this problem. Our article is organized as follows. 
In Section 2 we build up an effective motion equation of disorder heteropolymer to describe the different 
monomer alternation along the heteropolymer chain | p8| . The form of effective Hamiltonian of the problem 
under consideration is obtained in Section 3 to avoid discrepancies of different methods. Section 4 is devoted to 
consideration of a SUSY scheme ||l^ that is optimal for study of the disorder heteropolymer in weak segregation 
limit po| . In Section 5 a SUSY correlator is considered to examine memory and non-ergodicity effects (Section 
6). Finally, Section 7 is devoted to discussion of obtained results. 



2 Effective motion equation 

Usually, the SUSY method is based on the dynamic Langevin-type motion equation [p^ . But covalent bonds 
in the polymer chain make such way inefficient because dynamic theory in polymers is much more complicated, 
than statistical mechanics of usual many-body systems ||2^. Therefore, it is essential to obtain the effective 
motion equation instead of dynamic one. 

Let us initiate with directed polymer that is represented by the Gaussian chain with probability ^'(R, N) to 
find the end-point N at coordinate R. As is known the function ^'(R, N) obeys the Schrodinger-type equation 



with imagine time — iiV |21 



d'i'/dN ^ {Dd^/dli^ -UiK^N))^, (1) 

where number N ^ 1, D = is effective diffusion coefficient determined by the Kuhn segment length b, 
U (R, N) is external field. In the limit oo a solution of Eq.(l) can be represented in the form of functional 

integral over dependence r(n) of chain coordinate on number of internal monomers: 

^'(R,Ar) = J eM-SiiN{rin)}/2D)drin). (2) 

Here conventional action 5(R, N) = S'RAr{r(n)} — Lo{r{n))dn corresponds to the fixed end-points r(0) = 0, 
and r{N) = R is determined by Lagrangian of the Euclidian field theory jl^, |2^ 

where the effective kinetic energy, for which the continuum limit r{n + 1) — r(n) dr(n)/dn is taken into 
account, presents the covalent bonds between monomers of polymer chain |2l|| . Inserting of ansatz (|^) into 
Eq.(|l|) gives the Jacobi-type equation 

dS ^d'^S 1 fdS^^ 



After introducing generalized momentum p = dS/dH. and total derivative dp/dN = dp/dN + {pd/dTVjp, 
Eq.(4) takes the linear form of Burgers equation: 

dp/dN = D (a2p/9R2 + 2dU/dR) . (5) 

The above pointed relation of Eqs.(0)-(||) is a well-known fact in the theory of directed polymers, stochastic 
growth and kinetic roughening phenomena (see |2^]). 

The basic observation for our purpose is that the Schrodinger-type equation (|l|) takes the form of the 
Fokker-Planck equation ||2^ 

2 



dP d 



for distribution probability 



P(R, N) = *(R, N) exp{-F(R)/2D}, (7) 
which dependence on R is determined by effective potential 
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y = - j FdR. (8) 

Corresponding force F relates with the initial potential U in Eq.(^ as follows: 

1 , 1 9F 

As is well known ]2^ , Eq.(^) determines the probability to realize solution of corresponding Langevin-type 
equation 

5R/a7V = F(R, TV) + C(A^) (10) 
for stochastic variable R — R(A^). Here the separated stochastic force C, is fixed by the white-noise conditions 

(C(iV))=0, {aN)aN')) = 2D5{N -N'), (11) 

where the angular brackets denote average with respect to the distribution (^ . 

In going from the above considered case of directed polymer to the main object of our interest, the random 
heteropolymer AB, the coordinate R of end-point N turns into a stochastic Ising variable 6'(n), where 9{n) = 1 
if n-th segment is of type A and 9{n) = — 1 otherwise. Corresponding to the quenched disorder in fixed sequence 
of different type segments, the law 0{n) of monomer alternation along the chain is described by master equation 
type of that determines the Glauber dynamics [ p5[ . Corresponding sequence correlator a{n)a{n') for effective 
spin <j{ri) = 9{n) — 9{n), being deviation of microscopic value 9[n) from average 9{n), takes the form lilt] 



a(n)a{n') = C2 exp(— |n — n\/l), 

C2^4/(l-/), /^(1/2)(1 + 0R), (12) 

where overbar denotes the averaging over composition (quenched disorder). Here two sequence characteristics 
appear being correlation length I and fraction / of type-A monomers. 

Stochastic variable cr(n) that possesses of the correlator (^2|) with exponential form, is governed by the 
effective motion equation 

dcr/dn = -a II -t- s(n), (13) 
where the stochastic source s(n) is the white noise: 

s(^=0, s{n)s{n') = 2C2r^S{n - n'). (14) 
Relation between the microscopic value cr(7i) and stochastic (5-correlated variable s(n) is given by equation 

a{n) = / e-("-™)/'s(m)dm, (15) 



that is the solution of Eq.(p^). On the contrary to the colored noise cr(n), the white noise s(n) possesses of the 
Gauss distribution function 



P{s{n)} = (AC^Tr/l)-'/' exp j"^ s\n)dn^ , 



(16) 



that determines quenched disorder with intensity 46*2 / ^. Respectively, the local averaged field 

77(r, n) = (4C2)-i/V(n),5(r - r(n)) (17) 

represents the order parameter (here and below the monomer volume takes to be equal unity) . 

Apart from usual terms type of that in Eq.(|lO|), effective motion equation for the field ( p^ ) must contain 
the inhomogeneity contribution Dd^rj/dv^ (as the first term in r.h.s. of Eq.(^)) that takes Fourier transform 
—Dk^rjy^. As a result, the motion equation for the Fourier transformation 



[n) = / r;(r, n)e-*""dr (18) 
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of the field (|T^) takes the Langevin form: 



aW5« = -(afc)V-5W/977^ + Ck. (19) 

Here, as above, the continuum approximation for effective time n ^ 1 is used, the characteristic distance a = 
= Q-^/'^h is determined by the Kuhn segment length 6, and determination of effective force /k = —dTi/drj^ 
(cf. Eq.(||)) are taken into account. The white noise Ck — Ck{n) is defined by conditions type of Eqs.(pl]): 

(Ck) = 0, {Q{n)Ck'{n')) = <5kk"5(n - n'), (20) 
where angular brackets denote average over thermal disorder. 



3 Effective Hamiltonian 

To obtain the effective Hamiltonian Ti-lr]} = fl'{m} ||2^ in Eq.(p^ the thermodynamic potential il'{m}, being 
the average with respect to both the conformation and the sequence scattering, should be determined as a 
function of the averaged order parameter 

m(r) ^ ^(rKr,n)) = (4^2)-!/^ ^ (a(n)<5(r - r(n))). (21) 

n n 

With this aim, let us write the partition function in the form of the functional integral [|l9l 



Z = y (5TO(r) exp |c2X J m^(r)drj> x 




r(n)) - l| (5 |^(4C2)-i/2CT(n),5(r - r(n)) - m{r)jj. (22) 

Here x > is the composition Flory parameter, the first (5-function takes into account incompressibility condi- 
tion, the second one reduces to determination (^) of order parameter m(r). Further, one follows to represent 
the (5-functions as functional Laplace expansions over auxiliary fields Jp, J„i. Then the average expression 
in Eq.(22) takes the exponential form with the exponent J{Jp + Jmm)dr — ^{Jp, Jm}, where the last term, 
being averaged over sequence and configuration sets, gets conventional free energy at given fields Jp, Jm- The 
steady-state magnitudes Jp, Jm of these fields are determined by conditions 6fl/6Jp = 0, Sil/SJm = —m. In- 
serting expression for Jp, Jm to functional f2{Jp, Jm} gets the thermodynamic potential defined by equation 
Z = j (^TOk exp(— ri'{mk}), where (see more detailed Refs.|, |l|,|27|) 

f7' = ^rk|mkP + i^i«kk'kknmkf + J v{r)dr, (23) 

k kk' 

Ti, = T+{akf, T = r^~C2X- 



Here the inhomogeneity contribution is taken into account in the quadratic term, averaging procedure with 
respect to the quenched disorder (|l^) results in appearance of Z~^-terms, double k-dependence of the kernel 
it'kk' = 4iV~^(Za)~^(k^ -t-k'^)^^ is stipulated by averaging over the distribution (||). It is worth mentioning that 
expression ( p3| ) agrees with Refs.]^, 0, whereas in Refs.[^, |2^, where replica approach was used, the second 
summand of Hamiltonian ( p3| ) with the opposite sign was obtained. To correspond the self-action effects the 
integrand in last term has the usual expansion form [|l], |ll| 



V = -{^I/3\)m^ + {X/4\)m'^■, 

fi = i2C3c:^^^^r\ A = 24(1 + 5Ci/C2)r^, (24) 

C2 = 4/(l-/), C3 = |l-2/|. 
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As mentioned above, this method is based on the Eqs.jjl^) that allows to express quenched disorder correlator 
within the transfer matrix approach. If one uses the more popular replica method, then the field Jm and the 
order parameter m must take a replica index a under which in the Hamiltonian ( p3| ) the summation from 1 to 
n ^ should be made Then the quadratic contribution takes the form 

A„„(k)|TO„(k)|2 + \Y1 ^a/3(k)m„(k)TO^(-k), (25) 

where in the limit n ^ for coincided replica indices is Aaa{}^ — > 2Tk. As it was found out in spin glasses, 
the peculiarity of the systems with quenched disorder is governed by the hierarchy of the phase space that 
is characterized by random overlapping of the different replicas Q. Therefore overlapping parameter ylQ^(k) 
in the second term (^5|) is the stochastic variable under which the average should be made. Let us take the 
corresponding distribution in the simplest Gauss form 
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(26) 



where dispersion Wkika = cr'^{la) "^N ^(k^ + k|) ^ is given by parameter a (see ^). Then after averaging of 
the partition function Z = J £)mk exp(— Xl'lmk}) the second term in Eq.(p5|) takes the form 

^ ^ X! "kik2"^a(kl)TO^(-ki)ma(-k2)TO^(k2). (27) 

klk2 

As a result the thermodynamic potential ( p3| ) in the replica form is 

^' = ^rk|m„(k)|2 + i ^ ?i;kik2l™a(ki)P|TOa(k2)P + 



2 

ka kik2 



J v{ma)dr- ^ ^ Mkik2W„(ki)m/3(-ki)ma(-k2)m/3(k2). 



(28) 



klk2 



The Eq.(^) removes the above mentioned contradiction in the choice of the effective Hamiltonian form: the 
positive contribution of the second term is governed by the intra-replica interaction (see [ pT| , Q ) , whereas negative 
contribution in [ pS] , |^ is determined by replica overlapping. The distinction of mentioned terms makes itself 
evident in the fact that the former leads to renormalization of the value Tk, whereas the latter corresponds to 
the memory and non-ergodic effects. 

In order to carry out such renormalization (see |2^) using of the mean- field approximation it needs to replace 
one of the multiplier |mQ,(k)p in the second term of Eq.(p8|) by the bare Green function Gko that corresponds 
to w = 0, Ukik2 — and is determined by equation 

=r + 2a'(fc-fco)^ (29) 

independent on replica number a. To determine the parameters r, ko one needs to substitute the Eq.(|29|) into 
the corresponding Dyson equation 

Gko =T-k + 5I^"kk'Gk'o. (30) 

k' 

Then after integration over wave vector k' one obtains 

r = r + (3/47r)/-2(2r)-i/2^ k^^ = 2Tr^/^l{2r)^/^a. (31) 
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According to the first equation the positive determined parameter r increases smoothly with growth of bare 
parameter t (r ~ at t < 0, |r| ^ 1 and r ^ t at r » 1 (see Fig. la)). It means that without self-action 
effect, random heteropolymer is stable with respect to the microphase separation The second equation ( |3l| ) 
means that space period A = 27r/fco depends on renormalized thermodynamic parameter r, as it is inherent in 
random copolymers According to Fig. lb with increasing of r the value A grows monotonously from to oo. 
In so doing the greater the correlation length I, the faster change of the period A near the point t = 0. 

The final form of the effective Hamiltonian of random heteropolymer follows from the thermodynamic 
potential (|2|) renormalized under acting of the fluctuations 

H = Yi\\r]a{k)f J w(?7a)dr-^ ^ Ukik2?7a(ki)?7/3("ki)??„(-k2)?7/3(k2), (32) 

kct a klk2 

T-k = r + 2a^{k - fco)^ Wkik2 = <jHla)-^N-\]4 + \4)-\ 
Here the kernel v{r]a) is determined by Eqs.(p^), where m is replaced by rja- 



4 Supersymmetric scheme 

To build a SUSY scheme on a basis of effective motion equation (|l^ ) let us introduce the generation functional 

in 



\ dn Si]^ 



dot 



TT = TrT + 2a'(^-fco)V, 



(33) 



averaged over noise Ckii^)i where J-function considers the motion equation (19), the determinant is Jacobian 
of transfer from ^k to rji^. Then, the functional Laplace representation is used for ^- function that introduces a 
ghost field ipi^{n). To attach exponential form for the determinant in Eq.(|33[) Grassmann conjugate fields ipk{n), 
V'k('T-) should be used, that obey the conditions |Q 



j Si) ^ j 54) ^ j ^p^pS^ip = 1, = S^/jS4), (34) 

where figure brackets denote anticommutator. Then, assuming that the averaging over Ck("-) in Eq.(|3^) is 
determined by Gauss distribution with variance 1 (see Eqs.(^o|)), the standard form of the partition function is 
obtained: 



Z{r]} = J P{-q, if), ipjdipS^ij), 

P{r],ip]^,ip} = e-Kp{-S{r],ip]^,ip}), S= Ldn, (35) 

JQ 

L= f [{ipri - U - (^72) + {H'{7^}v - i)n"{T^}i))]dv. 



Here point denotes derivative over "time" n, the prime denotes the functional derivation with respect to the 
field (p7|), the last expression is taken in r-representation. 

The last expression in Eqs.(p5|) takes the simplest form |Q 

L=f A{^)d^d, A = ^(P$^)(2?$k)+W{$k}, d^i} = dM^ (36) 

k 

if one introduces SUSY generators 



and SUSY field 
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$ = ?7 + V'l? + iJV' + ^^V. (38) 

where Grassmann coordinates -d, d obey the same relations ( p^ ) as for the fields -0, "0- Here functional Ti-i^} has 
the same form as the effective Hamiltonian (p^), where order parameter mta is replaced by supcrficld <i>k(i?)7 
Eq.(|38|). In this case replica index a is removed by Grassmann variable § that is the formal reason of the replica 
and supersymmetry methods identity. The advantage of the latter method is in the Grassmann conditions (jsj) 
that correspond to the limit n ^ in the replica method. 

According to consideration |2^, the physical meaning of the components of the SUSY field ( |38|) is 
as follow: tp is the most probable value of fluctuations of the field conjugated to the order parameter 77, and 
combination 0-;/' determines density of sharp interphases. So, using the 4-component SUSY field ( ^8|) corresponds 
to the strong segregation limit ||^, |^ . Further we shall consider more simple case of the weak segregation limit 
where t/'V' = 0. Then, the SUSY field (|38| ) is reduced to the 2-component form 

(t) = T] + 9ip, (39) 
where self-conjugate nilpotent variable 9 = is introduced. Respectively, Lagrangian (pq) takes the form 



L= f A{(b)d9, A = Y,KD^k + n{(^k} 

k 



The motion equation for the nilpotent field (^^ that corresponds to Lagrangian (^) reads 

D<j>^ = -m/Scj^l (41) 

In the component form it leads to the equations for the order parameter 77(71) and the amplitude of the more 
probable fiuctuation (p{n) of the conjugate field (see pC[|). 

5 Correlators 

Now let us consider SUSY correlator 

Ck(7i, 9- n', 9') EE (0^(71, 0)0k(n', 9')) . (42) 
Multiplying Eq.(^l|) by value 0^ and averaging within zeroth approximation (v = 7i = in Eqs.(|3^)) one gets 

B 

'k -r 

Here conventional frequency v denotes Fourier transformation over "time" being the monomer number n, for 
example 

a = / C{n)e'''"dn. (44) 







The expression (|4^ ) has characteristic structure with respect to combination of the nilpotent variables 9, 9', 
that is inherent in, obviously, not only zeroth approximation but arbitrary supercorrelator. In this connection 
it is convenient to introduce basis supervectors 

A{9,9')^9, B{9,9') = 0', T{0,9')^1, (45) 
that have functional production 

X{9,9')= [ Y{9,9")Z{9",9)d9" (46) 
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for any vectors X, Y, Z. It is easy to see, that the basis SUSY vectors ( [45[ ) obey the following multiplication 
rules: = A, = B, BT = T, TA = T, other products are zero. Because the set of vectors A, B, T is 
complete, it is convenient to expand any SUSY correlator over this basis: 

C = G-A + G+B + ST. (47) 

Here and below subscripts k, v are suppressed for brevity. Using Eqs.(^9|), (^2|), one obtains for coefficients of 
the expansion 

G_ = (w*), G+ = {rj*ip), S={\rj\'). (48) 

So, G± represent advanced and retarded Green functions and S is the structure factor. In accordance with 
Eqs.(^), (^5|), (^7|), these functions within zcroth approximation take the form 

G^i^ = ir±i,y)-\ 5(°)=GV°^GL"^ = (r2+^.2)-i. (49) 
The Dyson equation for SUSY correlator ( ^ ) is as follows |^9[ ^ 

= (C(°))-i -uC-S. (50) 

Here S is self-energy function, u = a'^ {2N)^^ (lako)^'^ is the typical value of the interreplica overlapping Ukk' at 
k = k' = kg . By analogy with the expansion ( p7| ) the self-energy superfunction S that describes the self-action 
effects is 

S = S_A + S+B + ST. (51) 
Then using Eqs.(p9|), the Dyson equations for component G±, S take the form 

-f uG± = (r±ii/) - E±, (52a) 

S^(l + 27rC2r'^S{iy) + I])G+G_(1 - uG+G-)-\ (52b) 

Here (5-term is caused by quenched disorder. 

To complete the system (52) it needs to express the components of the self-energy superfunction by super 
correlators. Using the SUSY perturbation theory with accounting of the cubic and quartic anharmonicities (|2~ 
for matrix element of the self-energy function one obtains 



S(z, A = ^ (G(z, z'))' + ^ (C(z, z')f . (53) 

Here z = {r,n,6} denotes the set of coordinate r, "time" n, and nilpotent coordinate 6. It is important that 
the multiplication rule in Eq.(^3|) is not the same as for the functional product (^). Here the ordinary product 
should be used according to the following multiphcation rules ^ |2^ AT ^TA = A, BT = TB = B, = T, 
other products are zero. Then coefficients of the SUSY expansion (|3l[) take the form 

S±(n) = (^Ai' + Y^W) Sin)G±in), (54a) 

E(n) = i(^^2 + ^5(n))52(n), (54b) 
where r-representation for macroscopically homogeneous system is used. 
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6 Memory and non-ergodicity effects 



Following to Edwards and Anderson pC| ] let us introduce the composition memory parameter q = {rj{n = 
N)ri[n = 0)), which value determines the long-range correlation in the different monomer alternation along the 
whole heteropolymer chain. Moreover, the non-ergodicity parameter A = go — <? that is difference between the 
isothermal susceptibility go = G_ {v — 0) and thermodynamic value g = G_ {ly 0) is used. Then, the main 
correlators acquire the elongated form: 

G±H = A + G±oM, S{n) = q + So{n), (55) 

where index denotes the components corresponding to ergodic system without memory. Substitution of 
Eqs.(ll) into Eqs.(54) gives 

S±(n) = (^fi^ + y <7^ q (A + G±oH) + S±o(n), 

E±o(n) = + A2q)5o(n)G±o(n) + ^S^oi^)G±o{n); (56a) 



E(n) 



A2 



(56b) 



Here the terms, being non- linear with respect to correlators G±o, Sq, are included into summands Sj-q, Sq, the 
terms that contain production S'qA « have been dropped and, finally, one separates out the first summands 
that disappear when parameter q goes to zero. It is characteristic that to determine the self-energy functions 
(56) the "time" -representation was used, whereas the Dyson equations (52) require "frequency" Fourier trans- 
formations of non-linear expressions (56). To avoid this difficulty the fluctuation-dissipation theorem is used 



5o(n^O) = G±o(i^^O) =5, 

^±o{^ ^ 0) - So(n ^ 0) = ^(^2 + X\)g^ + 



(57a) 
(57b) 



where Eq.(|56q) is taken into account in the last equation. 



Inserting Fourier transformations of Eqs.(p5[), (56) to Dyson equation (52b), within the z^-representation one 
obtains 



90 



2 1 ^ 2 



So 



(l + Eo)G+G_ 



l-[u+{fl^ + X^q/2)q]G+G- 



(58) 
(58a) 



The first of these equations corresponds to S-teim being caused by memory effects, the second one - to "fre- 
quency" 0. At ^ the characteristic production is G+G- g^, so that the pole of structure factor 



(58a 



2 ^ \ —2 

fJ- + y9 ) 9 = .9 



(59) 



deter mines the point of ergodicity breaking. By analogy, substituting of Eqs.(56a), (57b) into the Dyson equation 



( 52a ) and taking into account g = G-{v 0) for microscopic susceptibility and memory parameter we have 

1 - rg + ug^ + ^.g [(.g + qf ~ q'] + ^g [{g + qf - q^] = 0. (60) 
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The thermodynamic behaviour of random heteropolymer in the vicinity of the ergodicity breaking point is 
described by the system of equations (|5^)-(6C). By analogy with spin glass Eqs.(58), ( |60| ) play role of 
Sherrington-Kirkpatrick equations, and Eq.(59) determines the point of de Almeida-Thouless instability. On 
further analyzing of these equations one should distinguish macro- and microscopic values qo, q oi the memory 
parameter and corresponding susceptibilities go, g. The peculiarity of such a hierarchy is that microscopic 
values, which conform to the limit — > 0, are the usual thermodynamic parameters and depend on temperature 
(Flory parameter x)- The macroscopic values qo, go conform to the point ly — and only depend on quenched 
disorder parameter I. In the non-crgodic area macroscopic quantities take the values at the point of ergodicity 
breaking. 

Because the system of three equations (|58|)--(|60|) is insufficient for determination of four quantities qo, q, go, 
g, it must be completed by equation 



1 



ug 



1 



X 



C2l-'g' 



(61) 



that is microscopic analogue of the equation ( |58| ) obtained as a result of separating out of the singular (5-terms 
for the structure factor in the Dyson equation. As is known from the spin glass theory 0|, the hierarchy of 
such singularities, which correspond to a set of infinitely decreasing "frequencies" — > 0, is a key point of the 
non-ergodic systems. In our case, the equations ([5^), ( |6l| ) correspond to the point 1^ = and a minimal of these 
"frequencies" respectively. In turn, at = the equation reads: 



A*^ + ) 90 = 5o ^ 



(62) 



Eqs.(pq), (f30|)-(62) are the complete system for determination of the quantities qo, q, go, g- 



7 Discussion 

According to Eqs.(p8|), (p2) the macroscopic memory parameter is given by cubic equation 



{fj.y2 + X\omql^C2r- (63) 
The characteristic form of the dependence of the value qo on intensity of quenched disorder I are depicted in 



Fig. 2. For copolymers close to symmetric composition / = 0.5 (C3 <C C2) the first term in brackets of Eq.(63) 
is negligible, and the dependence qo oc /^/-^ takes place. In the opposite case of dilute copolymer, where / ^ 1 
(C2 <SC C3), one obtains the dependence go (x f l^^^. 

Simultaneous solution of Eqs.(^, ( |60| ) and ( |6^ gives the point of ergodicity breaking xo, which dependence 
on correlation length I is depicted in Fig. 3 (solid curve). A characteristic of this phenomenon is that non-zeroth 
value of xo appears above a critical value of correlation length I, and with further growth of I dependence reaches 
its maximum and then monotonously falls down. In so doing, ergodic region is located under curve Xo(0 ^^'^ 



contracts with correlation length growth. Condition dg/dx = —00 added by Eqs.(|60|), (61) gives the equation 

u + n\gf + q) + —{gf+qf^gf, (64) 

that defines Xf value of Flory parameter in the freezing point, under which the microscopic susceptibility g 
takes zeroth value (see Fig. 4). Corresponding dependence x/(0 on the correlation length / is depicted in Fig. 3 
(thin curves). A characteristic is that the dependence %/(/) is under ergodicity breaking curve Xo(0 ^^'^ ^^^^ the 
same form. The influence of the composition on the values of xo and Xf is shown in Fig. 3a. As we recede from 
the composition / = 0.5, the growth of the above mentioned parameters is observed. The more complicated 
behaviour is realized at the growth of the interreplica overlapping a (see Fig. 3b). At small magnitudes of 
the correlation length I, the growth of both xo and Xf values is observed with cr increasing, whereas at large 
magnitudes I, parameters xo and Xf decrease. 

Dependencies of the macroscopic 170 and microscopic g susceptibilities on the parameter x are depicted in 
Fig. 4. Under point of ergodicity breaking xo, these susceptibilities and corresponding memory parameters q, qo 
coincide. The dependence g{x) has a cut-off at the point X = Xf (below it, susceptibility g takes zeroth value 
that corresponds to the freezing state). Above the ergodicity breaking point xo: the macroscopic susceptibility 
go is constant and the microscopic one has smooth decrease (latter can be obtained by simultaneous solution 
of Eqs.(pG), (pH)). According to Fig. 4a, with moving away from the composition / = 0.5, the values of the 
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susceptibilities go, g in the freezing point and in the point of ergodicity breaking decrease but corresponding 
values of parameters xo and Xf increase. The dependence on the correlation length I is shown in Fig. 4b: with 
increasing of I ergodic area contracts as it should be. The influence of the interreplica overlapping is shown 
in Fig. 4c. With increasing of the corresponding parameter a, the values of the susceptibilities g{x) and (?o(x) 
decrease and hence the interreplica overlapping prevents the heteropolymer freezing. 

The influence of the thermodynamic parameter x on the microscopic memory parameter q is depicted in 
Fig. 5. The absence of memory below Xf is the characteristic feature of the freezing region. Non-zeroth value 
of q appears in the freezing point Xf and with further growth of the parameter x, the memory parameter 
monotonously increases. The step-like behaviour of the parameter q is inherent in the first order transition. 
Evidently, the physical reason of the mentioned behaviour is the fluctuation contribution into the thermodynamic 
potential of the heteropolymer. According to Fig. 5a, with moving away from the composition / = 0.5, the 
dependence q{x) becomes more slight. On the contrary, the growing of the correlation length / results in more 
abrupt growth of the memory parameter (see Fig. 5b). At last, one can see from Fig. 5c that influence of the 
interreplica overlapping cr above and below of the ergodicity breaking point occurs to be opposite. 

According to the Fig. 6, the non-ergodicity parameter A(x) grows monotonously with increasing x from the 
ergodicity breaking point xo- The deflection from the composition / = 0.5, the decrease of the correlation 
length I and the growth of the interreplica overlapping a result in weakening of the non-ergodicity effects (see 
Fig.6). 

In the search for a new polymers with predetermined features, the phase diagram plays the basic role that 
defines possible thermodynamic states at different values of the Flory parameter x and composition /. According 
to Fig. 7, such diagram has a concave form for both the freezing point Xfif) and the point of ergodicity breaking 
Xo(/)- A region of the large x adjacent to the composition / = 0.5 corresponds to the non-ergodic unfreezing 
state. With decreasing of the values x and |/ — 1/2 1, at first the system goes to the ergodic state and then it 
freezes. The comparison of the Fig. 7a and Fig. 7b shows that the increasing of the correlation length I results in 
expansion of the unfreezing and non-ergodic phases. On the contrary, from Fig. 7a and Fig. 7c one can see that 
interreplica overlapping leads to its contraction. 
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Captions 



Fig.l Dependence of the parameter r (a) and the period A (b) on r at different values of the correlation length 
I (curves 1, 2, 3 correspond to I = 0.5; 1; 10). 

Fig. 2 Dependence of the macroscopic memory parameter qq on the value of the correlation length I (curves 1, 
2, 3 correspond to / = 0.5; 0.3; 0.1). 

Fig.3 Dependence of the characteristic values of the Flory parameter in the point of ergodicity breaking xo 

(thick curve) and in the freezing point Xf (thin curve) on the correlation length /: a) at = and different 
values of the composition / (curves 1, 2 correspond to / = 0.5; 0.3); b) at / = 0.5 and different values of 
the interreplica overlapping a (curves 1, 2 correspond to cr = 0; 2). 

Fig. 4 Dependence of the microscopic g and macroscopic go susceptibilities on the parameter x- a) at cr = 0, 
I = 0.1 and different values of the composition / (curves 1, 2 correspond to / = 0.5; 0.3); b) at / = 0.5, 
cr = and different values of the correlation length I (curves 1, 2, 3 correspond to / = 0.05; 0.1; 0.2); c) at 
/ = 0.5, I = 0.1 and different values of the interreplica overlapping parameter a (curves 1, 2 correspond 
to (7 = 0;2). 

Fig. 5 Dependence of the microscopic memory parameter q on the Flory parameter x- a) nt a = 0, I = 0.1 and 
different values of the composition / (curves 1, 2 correspond to / = 0.5; 0.3); b) at / = 0.5, cr = and 
different values of the correlation length I (curves 1, 2, 3 correspond to I = 0.05; 0.1; 0.2); c) at / = 0.5, 
I = 0.1 and different values of the interreplica overlapping parameter a (curves 1, 2 correspond to cr = 0; 2). 

Fig.6 Dependence of the non-ergodicity parameter A on the parameter %: a) at cr = 0, Z = 0.1 and different 

values of the composition / (curves 1, 2 correspond to / = 0.5; 0.3); b) at / = 0.5, cr = and different 
values of the correlation length I (curves 1, 2, 3 correspond to I = 0.05; 0.1; 0.2); c) at / = 0.5, I = 0.1 and 
different values of the parameter a (curves 1, 2 correspond to cr = 0; 2). 

Fig. 7 Phase diagram of the disordered heteropolymer: a) at cr = 0, a = 1; / = 0.1; b) at cr = 0, a = 1; 
Z = 5;c)atcr = 2,a = l;Z = 0.1. Thick curve corresponds to the ergodicity breaking point and thin 
curve corresponds to freezing point; FE^ NE, NN correspond to freezing ergodic, non-freezing ergodic, 
non-freezing non-ergodic phases. 
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